Five models for myosin V 
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Myosin V was the first discovered processive motor from tiie myosin family. It lias therefore been 
subject of a number of mechanical, kinetic, optical and structural studies and now belongs to the best 
characterised motor proteins. This effort has been accompanied by a number of different theoretical 
models. In this article we give an overview of them and discuss what they have in common and 
where the open questions are. The latter include the existence of sub-steps, the process that limits 
the run length, the nature of backward steps, the flexibility of the lever arm and the state of the 
lead head. 
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I. INTRODUCTION 

Myosin V is a two-headed processive motor protein 
from the myosin superfamily [1], involved in differ- 
ent forms of intracellular transport |2l E]. In 1999 it 
was shown to be the first processive motor from the 
myosin superfamily [1]. Processivity means that a sin- 
gle molecule can move over a long distance (typically 
a micron) along an actin filament without dissociating. 
This discovery came as a surprise, because processivity 
was only known for microtubule-based motors (kinesin, 
ncd and cytoplasmic dynein) before. The second fea- 
ture that makes myosin V unconventional is its long step 
size around 35 nm |1]. After myosin V, other processive 
myosins have been discovered (VI [5], Vila [H], XI [7], 
X [8^ and according to some evidence IXb [9]), but the 
amount of knowledge gained on myosin V has made it one 
of the best, if not the best understood motor protein. 

The experiments have characterised it mechanically 
(usually using optical tweezers) [4l HO l [TT 1 [T2 l [13 1 [14 1 [TS l 
[T5], biochemically [T71|TSl|Til2U], optically HH |11I23] 
and structurally [21 HSj HI H?]. These studies have 
shown that myosin V walks along actin filaments in a 
hand-over-hand fashion [23l [28] with an average step size 
of about 35 nm, roughly corresponding to the periodic- 
ity of actin filaments jUnjUH HH JH], a stall force of 
around 2.5 pN [TT] and a run length of a few microns 
[TTl|29l|30]. Under physiological conditions, ADP release 
has been identified as the time limiting step in the duty 
cycle [m |T^ , which means that the motor spends most 
of its time waiting for ADP release from the trail head. 

Myosin V contains a head domain, which has a high de- 
gree of similarity with other myosins like myosin II (mus- 
cle myosin). However, there are important differences in 
the kinetics of the ATP hydrolysis cycle - notably the 
slow release of ADP, which is necessary for maintaining 
processivity. The head is connected with a lever arm (also 
called neck) , consisting of an alpha helix, surrounded by 
six light chains. This is followed by the tail domain, 
which is responsible for the dimerisation of the molecule 
and for cargo binding. 

Meanwhile the literature on myosin V counts about a 
dozen of reviews, concentrating on different aspects of 
myosin V, such as its cellular function and phylogeny 
[HIIIIIISIIMIIMIIMIES], regulation [Mj I3S], structure 
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FIG. 1: (A) Myosin V steps in a hand-over-hand fashion along 
actin filaments. (B) The most frequent step size is 13 actin 
subunits, but it can also be 11 and 15 subunits occur 



pr and stepping mechanism [3511351121110111111121113] ■ 
It is not the purpose of this article to review all important 
findings on myosin V in general, but rather to concentrate 
on theoretical models and those experimental results they 
aim at explaining. More specifically, we will critically 
review 5 existing theoretical models gH US HHl HZl HE] 
and discuss which of them is in agreement with which of 
the experimental results and where the discrepancies are. 
The structure of this paper is as follows. In Section jIll 
we summarise those experimental results that are most 
relevant for theoretical models. Section |TTT] describes 
some of the most important open questions in the field - 
particularly those that need cooperation between theory 
and experiment to be resolved. Section |IV| then reviews 
five theoretical models for the stepping kinetics and me- 
chanics of myosin V. Finally, we give a brief summary 
and try to project how the different approaches might 
converge in near future. 



II. MAIN FINDINGS 

The main findings gained with optical tweezers, fiu- 
orescent microscopy and chemical kinetics can be sum- 
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FIG. 2: Force- velocity relation of myosin V for different ATP 
concentrations. Data points are taken from Refs. [15II16II30] . 
The lines are guide to the eye. The relationship shows a rather 
fiat region between — lOpN and IpN. In that region, the ve- 
locity as a function of ATP concentration follows a Michaelis- 
Menten-like relation. This region is followed by an abrupt 
drop between IpN and 2pN. Above the stall force, the mo- 
tor is pulled backwards with a velocity that increases roughly 
linearly with load and is independent of the ATP and ADP 
concentration. 



marised as follows: 

1. Myosin V walks in a hand- over-hand fashion along 
actin filaments IJSl ^E\, as shown in Fig. [TlA. The 
lead head performs free Brownian motion around 
the common hinge while searching for the new bind- 
ing site [49, 50J. 

2. The step size roughly corresponds to the actin pe- 
riodicity of 35nm [5T] . The most precise way to de- 
termine the average step size is by measuring the 
pitch of the slightly helical motion of the stepping 
myosin V molecule around the actin filament 2T]. 
The distribution of step sizes can be determined 
from EM images [21] • The study shows the highest 
probability for 13 subunit steps (36nm), followed 
by 15 subunit and 11 subunit steps (Fig. flp). 

3. The force velocity relation, shown in Fig. [2] reveals 
the following phases: 

• Under saturating ATP concentrations (in 
the millimolar range), the velocity is about 
400nm/s and is nearly constant up to loads 
of about IpN [HEI]. 

For higher loads, the velocity drops rather 
sharply and fiattens near stall [HI [51] . 
Above the stall force, when the motor is 
pulled backwards by the load, its velocity in- 
creases roughly linearly with the applied load; 
it reaches around — 150nm/s [T5| to — 250nm/s 
[l4] under 5pN load and around — 700nm/s 
under lOpN [HJ. 
Forward loads only lead to a slight increase in 



Under low ATP concentrations, the veloc- 
ity becomes proportional to the ATP concen- 
tration and is independent of load between 
-lOpN and 2pN [13 (TSJ [5T] . For super-stall 
loads, the velocity becomes independent of the 
ATP concentration. 

The stall force of myosin V has been reported 
as 3pN [5, 2.5-3pN [TH], up to 2.5pN [55], 
but also lower values around 1.7pN [TJj. Be- 
cause of short run lengths, some care has to 
be taken when determining the force-velocity 
relation close to stall force, as pointed out by 
Kolomeisky et al. for kinesin |53j . 

Added ADP has a similar effect on the force- 
velocity curve as a reduced ATP concentra- 
tion [m [TBI I3D] • Added inorganic phosphate 
in high concentrations (40mM) causes a less 
significant, but detectable velocity decrease 

MM- 



velocity ^^. 



The second important source of data is the degree 
of processivity, usually measured with the average 
run length, i.e., the average distance a single myosin 
V motor moves along an actin filament before dis- 
sociating from it. There is still some discrepancy 
with regard to the run length of unloaded motors. 
The first measurement was carried out by Mehta et 
al. [1|, who estimated it as 2^m. Veigel et al. [T?| 
obtained 2.4 /xm, Ali et al. [51] 2.1 /xm, Baker et al. 
[50] measure values around 1 /im, while the data of 
Clemen et al. [14] are closer to 300nm. A possible 
explanation for this apparent discrepancy could lie 
in the fact that actin filaments were attached to the 
glass surface in the latter experiment, which pos- 
sibly hindered the motor in its natural helical mo- 
tion around the filament and caused a premature 
dissociation. Therefore, one can accept run length 
values around 2 /im for unconstrained motion along 
actin filaments. 

• The mean run length over a wide range of nu- 
cleotide concentrations, but without an ap- 
plied load, has been measured by Baker et 
al. [SO]- It shows a sHght decrease with a 
growing ATP concentration. Added ADP de- 
creases the run length, but only by about a 
factor of 2, which is reached at concentrations 
around 200 /iM. Under higher ADP concen- 
trations, the run length remains constant at 
around 400nm. Added 40 mM phosphate led 
to a reduced run length [30^ at zero load and to 
a lower load at which the molecule dissociates 

• Knowledge on the load-dependence of the run 
length is still rather incomplete. Clemen et al. 
[14] measure a largely constant run length of 
about 300nm under loads ranging between - 
5pN and 1.5pN. Under a super-stall load of 
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FIG. 3: Duty cycle of a dimeric myosin V motor. A head with ADP undergoes the power-stroke, then the second head 
containing ADP.Pi binds in the lead position, it releases Pi, then the rear head releases ADP, binds a new ATP molecule and 
detaches from actin. After that the cycle repeats with exchanged roles of the two heads and with the motor one step further. 



5pN, the average backward run length was 
SOnm. 

• The run length decreases with an increasing 
ionic strength ^U\ |5S] . It could be increased 
by adding (through mutation) additional pos- 
itive charge to loop 2, which increases the 
electrostatic attraction between the myosin V 
head and actin [55]. 

These mechanical data are directly complemented 
by kinetic measurements, which follow the transi- 
tions that take place when motor molecules with a 
certain nucleotide (ATP or ADP) are mixed with 
actin filaments and a different nucleotide in solu- 
tion. There are generally two ways of probing the 
myosin V kinetics - using single-headed (usually 
SI) molecules, or double-headed motors. Experi- 
ments on single-headed motors reveal the kinetics 
of the unstrained ATPase cycle. Experiments with 
double-headed motors, on the other hand, bring 
additional insight into the kinetics in states where 
both heads are bound and some transitions are sup- 
pressed (or possibly in some cases accelerated) by 
intra-molecular strain. 

Single- headed: 

Kinetic studies on single-headed myosin V 
molecules were pioneered by De La Cruz and 
coworkers. In Refs. [17l [20], a number of ki- 
netic rates have been determined. These include 
the ATP binding rate (both in free heads and on 
actin), the ATP hydrolysis rate, the phosphate re- 
lease rate on actin, the ADP release and re-binding 
rate (in free/bound heads), and the actin- myosin 
binding/unbinding rate in apo (state without any 
nucleotide) and ADP states. Trybus et al. ^S] re- 
port similar values, but with a large discrepancy 
with regard to the ADP release rate in free heads. 
Wang et al. [57], on the other hand, report a slower 
phosphate release rate. Forgacs et al. [58] measured 
the kinetics with deac-aminoATP. In summary, one 
can say that the kinetics of single-headed myosin V 
is well characterised, except for some short-lived 



states (such as the dissociation rate in the ATP 
state) or rare transitions (such as dissociation in 
the ADP.Pi state). 

Double-headed: 

Kinetic measurements on double-headed myosin V 
molecules in the presence of actin are crucial for a 
quantitative understanding processive motility and 
the coordination of chemical cycles. At the same 
time, they are much more difficult to interpret, be- 
cause one measures a superposition of contributions 
from the lead head, from the trail head and from 
molecules bound with a single head only. 

The first kinetic study on dimers was carried out 
by Rosenfeld and Sweeney [SS]. They found that 
myosin V with ADP, mixed with actin, binds with 
both heads and releases ADP from the trail head 
at a rate of 28 — 30 s~^ and from the lead head 
at about 0.3 - 0.4 s~^ Note that the ADP release 
rate from the trail head is somewhat (by about a 
factor of 2) accelerated in comparison with single- 
headed myosin V, while it is strongly suppressed 
on the lead head. This is attributed to the effect 
of intra-molecular strain, which affects the ATP re- 
lease both from the lead and from the trail head. 
If myosin V with ADP and Pi on it is mixed with 
actin, it releases phosphate from both heads at a 
measured rate of about 200 s^^. However, obtain- 
ing this rate involved some extrapolation, so the 
lower values obtained in later studies [60) 161] might 
be more reliable. But the important message of 
the paper is that the lead head releases phosphate 
quickly and then waits in the ADP state. ADP re- 
lease from the lead head was therefore identified as 
the key gating step in the chemical cycle. 

As a result of these kinetic studies, the duty cycle 
shown in Fig. [3] has emerged with a high level of 
consensus. 




No telemark state 



Telemark state 



FIG. 4: In the model without a telemark state (left), the lead 
head stays in the pre-powerstroke state before the trail head 
detaches. In the telemark state (right), both heads are in the 
post-powerstroke state, but the leading lever arm has a strong 
kink - resembling a telemark skier. 
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III. OPEN QUESTIONS 

Yet some of the important questions on myosin V have 
only been partially answered so far. In the following, we 
will discuss five major open questions. They are, on one 
side, important for setting up theoretical models. At 
the same time, most of these questions are quantitative 
in nature and most of them concern properties that are 
not directly observable in an experiment. Therefore, they 
can only be tackled efficiently by combining experimental 
results with theoretical models. 



A. Head-head coordination 

In order to maintain its processivity, a dimeric motor 
has to prevent simultaneous detachment of both heads 
from its track. This means that the chemical cycles of 
both heads have to be coordinated in order to stay out 
of phase. There is a high level of agreement that this 
coordination is achieved through intra-molecular strain. 
There is also undisputed evidence for lead head gating, 
which means that ADP release from the lead head waits 
until the trail head detaches from actin [59, 61, 62J. How- 
ever, there is also some evidence for much weaker gating 
in the trail head, meaning that the ADP release rate in 
the trail head is somewhat accelerated by the presence of 
the bound lead head IT2ll59l. 



Both effects are likely connected with the small (about 
5 nm for single-headed molecules) conformational change 
that takes place upon ADP release [TH [33] • Conversely, 
Oguchi et al. [SJ have recently directly measured the 
influence of an applied load on the ADP binding affinity. 
The strong deceleration of ADP release in the lead head is 
possibly also connected with the large main power stroke. 
The exact nature of this influence depends on whether 
the "telemark state" exists or not (see next section). 

In addition, coordination can take place during other 
steps, too. The elastic lever arm model shows that the 
attachment rate of the lead head is strongly suppressed 
before the trail head goes through a power-stroke j45| . 



FIG. 5: Sub-steps as observed by Uemura et al. (A) and Cap- 
pello et al. (B). Uemura et al. observe a substep of llnm, 
preceding the main step. Its size is independent of load, but 
its duration increases with load or upon addition of BDM. 
Cappello et al. observe one substep preceding the main step 
(only at low ATP concentrations) and one following the main 
step (independent of ATP concentration). 



B. Lead head state 

One of the most controversial questions on myosin V is 
whether the lead head is in the pre-powerstroke state, or 
in the post-powerstroke state with a strongly bent lever 
arm (also called "telemark" state, because the shape of 
the molecule resembles a telemark skier with the front 
knee bent), as shown in Fig. l4J Evidence for the tele- 
mark state flrst came from EM images [21] , although sub- 
sequent analysis revealed that molecules with a kinked 
leading lever arm are not as abundant as those with a 
straight one [2S|, and that the telemark state might not 
lie on the main pathway. Further evidence came from ex- 
periments with fluorescence polarisation |22[ 165] , which 
showed that about a third of labels attached on the lever 
arm did not show any detectable tilting motion. But later 
studies revealed that this was probably a consequence of 
the ambiguity in angle detection from one vector compo- 
nent alone and that all labels on the lever arm do indeed 
show tilting motion, as well as alternating steps [Ml [57] . 

However, it is still possible that the head domain of the 
lead head is in some intermediate state and that its lever 
arm is bent, but not as strongly as in the pure telemark 
state. 



C. Sub-steps 

An intriguing question is whether the motor moves 
its load in uniform steps, or do they contain short-lived 
substeps. Because of the noise in position determina- 
tion, these sub-steps are difficult to detect. The experi- 
ment of Uemura and co-workers [16] revealed substeps of 
ll-|-24nm, as shown in Fig. |5]A. Their length was inde- 
pendent of load, but their duration was increased both 
with an increasing load and upon addition of BDM (bu- 



tanedione monoxide). The duration was independent of 
the ADP and ATP concentration. 

The second important observation related to sub-steps 
comes from Cappello et al. [55]. This study employs 
travelling wave tracking, a technique that allows a sig- 
nificantly higher temporal resolution. It reveals two sub- 
steps of about 5nm: one immediately preceding and one 
immediately following the main step (Fig. pp). The du- 
ration of the substep before the main step was longer at 
a lower ATP concentration. Therefore, neither of the two 
sub-steps has a direct correspondence in the observation 
of Uemura et al. The substep following the main step 
can well be attributed to the diffusive search of the lead 
head for the next actin binding site. The explanation for 
the step preceding the main step is more difficult. The 
authors propose a hypothesis in which the substep rep- 
resents a reversible transition in the lead head that can 
take place after the trail head releases ADP and which is 
a necessary condition for ATP binding to the trail head. 
One question that remains open in this scenario is the 
mechanism by which ATP binding to the trail head could 
be triggered by a transition in the lead head. 

Sub-steps have also been observed in two somewhat 
different experiments. First, Dunn and Spudich |49j 
have labelled sites on the lever-arm with gold nanoparti- 
cles and tracked it with sub- millisecond time resolution. 
They found that the 74 nm steps were subdivided into 
49 nm+25 nm and that the intermediate state had a very 
high level of fluctuations, resulting from the free head 
diffusively searching for the next binding site. Depend- 
ing on the light chain presence, the lead head binding 
rate was either 60 s~^ or 180 s""'^. It was slowed down by 
the presence of BDM - different from the conclusion of 
Uemura et al. [TB] that the effect of BDM is a delayed 
power stroke. The lead head binding rate is somewhat 
slower (possibly because of the attached particle), but 
still of the same order of magnitude as the duration of 
the second sub-step (1 — 2 ms) measured by Cappello et 
al. |52j . A direct comparison of step sizes between the two 
experiments is difficult, because Cappello et al. tracked 
the cargo position, whereas Dunn and Spudich tracked a 
bead attached to the lever arm. 

Finally, as already mentioned, a sub-step of about 5 nm 
has been observed by Veigel et al. [T^[53] in single-headed 
myosin V molecules and identified as the conformational 
change upon ADP release. There are good reasons to 
believe that this sub-step is not directly observable in 
dimeric motors |45j . One argument is that ADP release 
shows a strong load dependence in single-headed motors 
[63LI64J . In dimeric motors, moderate backward loads up 
to '--^ 1 pN or forward loads show almost no effect on the 
stepping rate, which is (under saturating ADP concentra- 
tions) limited by ADP release IHHS]. However, Cappello 
et al. hypothesise that the conformational change upon 
ADP release in the trail head might be involved in trig- 
gering a reversible transition in the lead head, which then 
causes an observable sub-step [52, . 



D. How tight is the coupling between the 
hydrolysis cycle and the stepping? 

The Michaelis-Menten like dependence of the kinesin 
velocity on the ATP concentration [68] was seen as early 
evidence for a tight coupling between ATP hydrolysis 
and stepping (in fact, those experiments could not yet 
exclude additional futile ATP hydrolysis cycles). Similar 
experiments also revealed tight coupling between step- 
ping and ATP hydrolysis in myosin V under low loads 
[TU] . Therefore, early kinetic models assumed ATPase cy- 
cles tightly coupled to the stepping [111 [69] , which means 
that the motor always makes one step when it hydrolyses 
one ATP molecule and that it can only make a backward 
step upon re-synthesis of an ATP molecule from ADP 
and phosphate. 

Such reversible behaviour has indeed been demon- 
strated in the Fi- ATPase, which can either hydrolyse 
or synthesise ATP, depending on the direction of rota- 
tion |70| . However, experiments on both myosin V and 
kinesin have shown that tight coupling has its limits un- 
der high loads. For kinesin. Carter and Cross [71, found 
out that under high loads the motor is pulled backwards, 
for which it requires ATP binding. The observation that 
myosin V can be pulled backwards by a super-stall force 
was first made by Clemen and coworkers [14]. The back- 
ward and forward steps were further investigated in a 
subsequent paper [TS] , where they found out that forward 
steps steps involve ATP hydrolysis, while backward steps 
involve neither hydrolysis nor synthesis of ATP. A rather 
surprising consequence of this finding is that a high load 
(e.g., 5pN) can pull the motor faster when applied in the 
backward direction than in forward. 

So under sufficiently high loads, both myosin V and 
kinesin leave the tightly coupled cycle, but each in a dif- 
ferent way: kinesin can be pulled backwards in the pres- 
ence of ATP, whereas myosin V neither hydrolyses nor 
synthesises it. 



E. Source of elasticity 

In any mechano-chemical model for myosin V, there 
needs to be at least one elastic element, such as the lever 
arm. One reason is that some degree of ficxibility is nec- 
essary in order for the motor to be able to bind to the 
actin filament with both heads simultaneously. Further- 
more, many models propose that when a head commits 
its power stroke, some of the free energy released is first 
stored into elastic energy and later uses it to perform 
work on the load. Besides the lever arm, which could 
be an obvious elastic element, further candidates for the 
source of compliance lie in or around the converter do- 
main. In particular, the converter domain was shown 
to be the main source of compliance in myosin II |72j . 
However, the same study showed that the converter do- 
main was significantly stiffened through a single muta- 
tion, which means that the result is not directly trans- 



ferable between different myosin classes. 

Measuring the stiffness directly requires careful elimi- 
nation of other sources of compliance, such as the actin 
filament and the tail domain. Veigel et al. measured 
the stiffness of single-headed myosin V and obtained the 
value 0.2pN/nm |63j . A comparison with smooth mus- 
cle myosin led to the suggestion that the elastic element 
is located close to the pivot point of the lever arm. Of 
course, this conclusion cannot be regarded as definite, 
because the degree of mechanical similarity between the 
two types of myosin is not known. 

On the other hand, the existence of a telemark state 
would suggest that the lever arm has a very compliant 
point (sufficiently compliant that it can be bent by the 
power-stroke), possibly around the 4th IQ domain. MD 
simulations |73j support the possibility of such a kink, 
but are not yet capable to precisely predict the bending 
energy connected with it. 



IV. THEORETICAL MODELS 
A. Overview 

Theoretical models can be employed to describe the 
dynamics of a motor protein on different levels. First, 
there are discrete stochastic models that describe the mo- 
tor with a collection of chemical states that correspond to 
different positions along the track |75j . The motor pro- 
tein is treated as a "black box" , the only observable being 
the tail- or the cargo position. This approach is partic- 
ularly useful for interpreting experimental results gained 
with optical tweezers, such as force- velocity relations un- 
der different ATP, ADP and phosphate concentrations. 
It was first applied to myosin V by Kolomeisky and Fisher 
[33j. The basic version of a discrete stochastic model al- 
lows no events beyond the regular chemical cycle, but it 
can be extended to allow parallel pathways [13 HH |7S] , 
as well as dissociation from the actin filament. Coordi- 
nated hand-over-hand motility is one of the assumptions 
of such models, rather than being its outcome. 

The next level of models uses a similar phenomeno- 
logical approach to describe each individual head. The 
properties of a dimeric motor are then derived from those 
of the two heads. This approach is similar to that taken 
by muscle myosin models [TTJ [751 175] , which deal with a 
large number of myosin heads, assembled into thick fila- 
ments. Such models combine the head domain, described 
with a few discrete chemical states, with an elastic ele- 
ment that connects it to the other head and to the tail. 
We will therefore use the term mechano- chemical m,odels 
for this category. Two models for myosin V that fol- 
low this approach are described in Refs. [l5] and [15] . 
In theory, the mechano-chemical models can be mapped 
onto corresponding discrete stochastic models. A state in 
the resulting discrete stochastic model is then described 
with a combination of the lead head state, the trail head 
state and the relative distance of the two heads on the 



track. Without further simplifications, the number of 
such states is therefore usually too high to be studied 
analytically. 

Finally, there are molecular models for motor proteins 
that attempt to understand the mechanism of energy 
transduction in microscopic detail. A full scale molec- 
ular dynamics simulation over millisecond time scales on 
which the chemical cycle takes place is still well beyond 
the available computational power, but it can be used to 
probe essential processes in the course of ATP hydrolysis 
[80l [8l] l82] . Also, several approximative methods have 
been employed. The dynamics of the power stroke has 
been investigated with molecular kinematics using the 
minimum energy path method (MEP) [S3J [M]. While 
these computational models represent an important con- 
tribution to our understanding of the microscopic mech- 
anism by which molecular motors work, they are cur- 
rently not yet precise enough to make predictions on free 
energy differences between different states or transition 
rates between them. Therefore, quantitative models still 
need to be "reverse-engineered" to a large extent from 
experimental data. 

There are other noteworthy models that deal with spe- 
cial aspect of myosin V, but will not be the subject of 
our review. Here we should briefly mention the work of 
Smith, which relates the duty ratio of myosin V to its 
processivity |85j . Further, there is a study by Schilstra 
and Martin [86] who investigate how viscous load can de- 
crease the randomness of the walking molecule. Finally, 
a number of models deal with the collective dynamics 
of motors in general [87l [88], and those findings could 
be relevant for the dynamics of myosin V transporting 
vesicles, too. 

In the following, we will give an overview over 3 dis- 
crete stochastic models and 2 mechanochemical ones. In 
particular, we will identify the common points and the 
differences between them. Table [ij shows an overview 
which model properties are in agreement with which ex- 
perimental results and where the discrepancies are. 



B. Discrete stochastic models 



1. Kolomeisky and Fisher 



Kolomeisky and Fisher [ll] (Fig. ^K) presented the 
first theoretical model describing myosin V. It uses 2 
states per step. The intermediate state corresponds to 
a substep of about 13nm, measured from the ATP bind- 
ing state. The model parameters were determined by fit- 
ting the calculated force-velocity relations to data from 
Ref . [4] . The conclusions of the study are that both the 
ATP binding rate is load-independent and the ADP re- 
lease rate is only weakly load dependent - which implies 
that most of the load dependence is contained in their 
reverse rates. At super-stall forces, the model shows 
backward steps, but they require ATP synthesis. Be- 
sides fitting the load velocity relations, the model also 



TABLE I: Comparison of five theoretical (kinetic and mechano-chemical) models. "+" denotes that the given experimental 
observation is in agreement with the model, "o" means partial or qualitative agreement, "— " denotes disagreement, "n/a" (not 
applicable) that the given aspect is not included in the model or that no analysis is available that would allow a comparison. 
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*Model shows backward stepping, but only through ATP synthesis 

'Velocity too high; dependence on ATP concentration not available 

■fin the refined model [74] no backward stepping takes place above 3pN load. 

S Velocity too low 

^Model shows a substep upon trail head release, but predicts an incidence that increases with the ATP concentration. 
** Histograms of position changes show intermediate peaks, but no analysis is available which would attribute them to any type of sub-steps. 
ttModel shows a substep upon trail head release, but it does not depend on ATP concentration 
1-1- Model process! ve only up to 3pN forward load 



makes predictions on the randomness parameter of the 
trajectory [5^, defined as 



,. (x^W) - {x{t)f 

r = lim , , .. 

t^oo a {x[t)) 



(1) 



"where a denotes the step size and (...) an ensemble av- 
erage. The value of r is 1 for a stepper "with a single rate 
limiting step. A higher value of indicates the presence of 
back"ward- or multiple steps, "while r < 1 means that sev- 
eral rate limiting transitions per step are involved. For 
myosin V, the randomness parameter has not yet been 
measured to date. 

Prediction of substeps of a similar size as later observed 
by Uemura et al. [16' is a particular merit of this model. 
However, it should be mentioned that the observed sub- 
steps had a duration around 10ms. In the model, the du- 
ration of substeps corresponds to the ADP release time, 
"which is around 80ms. This discrepancy is an indication 
that t"wo states per step are not sufficient if one "wants to 
describe the myosin V motion in this amount of detail. 



2. Skau, Hoyle and Turner 

Skau et al. [ITj (Fig. l6]3) developed a more elaborated 
myosin V model, "which still belongs to the category of 
discrete stochastic models, but includes t"wo parallel cy- 
cles - one that leads to a for"ward step and a futile one 
that does not. In addition, one dissociation path from 
actin is included. The model does not include back"ward 
stepping "without ATP synthesis. Substeps of 25-t-llnm 
are included in the model, but in reverse order compared 
"with the model by Kolomeisky and Fisher [H] and the 
experimental observation by Uemura et al. |16j . "where 
ATP binding is followed first by a short (12nm) and then 
the longer (24nm) substep, rather than the long step fol- 
lowed by the short one. 

The model fits well the force- velocity relations, al- 
though the velocity starts to drop around zero force al- 
ready, rather than being nearly fiat up to about IpN, as 
experimental data suggest. Because the model does not 
include backward stepping by other means than revers- 
ing the ATP hydrolysis, it fails to reproduce the velocity 
at super-stall loads. 

The calculated run length increases inversely propor- 
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FIG. 6: Four kinetic models with different levels of complexity. (A) Model by Kolomeisky and Fisher [J^ (B) Model by Skau, 
Hoyle and Turner [47]. The right part shows the states numbered 1-7. Note that states 1,3,4,5 and 6 have the same load 
position. (C) Model by Baker et al. [30] (D) Model by Wu, Gao and Karplus j48j. 



tional with a decreasing ATP concentration. It also in- 
creases with an increased ADP concentration - the latter 
is inconsistent with the measurement of Baker et al. 1301. 



3. Wu, Gao and Karplus 

Wu et al. [H] (Fig. JBt)) recently proposed a model 
that contains more kinetic details, but no external force 
(load). In that model, each head can be found in the 
empty state, in the ATP state, in the detached ATP 
state, in the weakly bound ADP state and in the strongly 
bound ATP state. The trail head can either be in the 
empty, ATP, or strongly bound ADP state. The lead 
head can be in the strongly or weakly bound ADP state. 



The chemical cycle can follow three different pathways 
(mainly depending on when the weakly- to strongly- 
bound transition in the leading head takes place) . Upon 
trail head detachment, the tail moves by 12nm, corre- 
sponding to the substep from Ref . [16] . This substep oc- 
curs only on one of the three kinetic pathways, which is 
in agreement with the experiment that showed substeps 
only for a fraction of steps. There are also two different 
termination pathways, which lead to the dissociation of 
the whole dimer from actin. 

The main results of the model are zero-load velocity 
and run length as a function of ATP and ADP con- 
centration. The calculated run length shows a decrease 
with the ADP concentration, which is in good qualita- 
tive agreement with Ref. [30j. The run length increases 
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somewhat with an increasing ATP concentration, which 
is different from the conclusion of Baker and coworkers 
[5U] . although, given the uncertainty of data points, not 
necessarily inconsistent with their measurements. The 
model shows that the velocity at low ATP concentrations 
is very sensitive to ADP, which could possibly provide a 
reconciliation of the large differences between different 
experimental results. 

A question that is beyond the scope of a purely kinetic 
model, but will have to be answered when mechanical as- 
pects are included, is how the intra-molecular strain can 
have such a strong influence on the weakly- to strongly- 
bound transition in the lead head. In the model, the tran- 
sition rate decreases almost 4-fold when the trail head 
releases ADP. Another puzzle brought up by the model 
is why the dissociation rate is 30 times higher in a state 
where both heads are bound to actin (one weakly and one 
strongly) than in a state where only one head is weakly 
bound to actin and the other one is free. While this seems 
inevitable for an explanation of run length data by Baker 
et al. |30j , it is not easy to think of a mechanical model 
in which myosin V, bound with both heads, is more likely 
to dissociate than when bound with a single head. Inter- 
nal strain might potentially provide an explanation, but 
only if the dissociation takes place in a single step and 
not with one head after the other. 



C. Mechano-chemical models 

The main difference between the models listed above 
and the mechano-chemical approach is that the latter 
uses a single myosin V head as the model unit and de- 
rives the properties of a dimer. This reduces the num- 
ber of parameters and model assumptions. The model is 
specified by the states of an individual head (typically 4 
or 5 bound states, plus a number of free states), by the 
mechanical conformation of the molecule in that state, by 
the transition kinetics and by the interaction between the 
two heads that is mediated by the mechanical properties 
of the lever arms joining them. In practice, however, even 
the kinetics of a single head contains a number rates that 
have not yet been measured directly and that need to be 
determined from dimer properties, such as force- velocity 
relation and processivity. So far, two myosin V models 
have been developed on these premises: the model by 
Lan and Sun [46, and the elastic lever arm model |45| . 
Although they share the basic principle, there are signif- 
icant differences in the implementation. 



1. Lan and Sun 

Lan and Sun [JB] proposed a model based on a cy- 
cle with 5 chemical states, all of which can in theory be 
bound to actin or free. The transition rates are taken 
from kinetic studies or guessed on the basis of other 
myosins where no data are available. Nonetheless, some 



assumed rate constants differ significantly from generally 
accepted values. For instance, the detachment rate of 
a head from actin in the empty state is assumed to be 
0.16 s-i, while the value from Ref. [T7] is 0.00036 s'^ 

The elasticity of the two lever arms is modelled as 
highly anisotropic. It combines an in-plane component 
and a very stiff azimuthal component. The in-plane 
elasticity corresponds to a persistence length of 120nm, 
or a spring constant of 0.05pN/nm, measured at the 
tip. This is consistent with an earlier estimate based 
on tropomyosin [90], but much lower than the stiffness 
measured with optical tweezers (0.2pN/nm [53] )• The 
azimuthal component is modelled in an ad hoc way and 
is assumed to be very stiff. If the model molecule binds 
to two adjacent sites on the actin filament, the bending 
energy is 40 fc^T, compared to 15 ksT in the model with 
isotropic lever arm elasticity [45 . 

In contrast with the model in Ref. [44 , the heads also 
possess internal fiexibility, meaning that the lever arm 
angle can fiuctuate to some extent while the head stays 
in the same chemical state. This allows the transitions 
to take place without the need for a simultaneous lever 
arm swing. 

The main result of the model is the force- velocity rela- 
tion for a fixed set of nucleotide concentrations. It repro- 
duces well the flat region at low loads and the abrupt drop 
around IpN, but there is a discrepancy above the stall 
force, where it predicts a backward speed well above ex- 
perimental values. No processivity calculations are made 
in the paper, so it is likely that the parameters will need 
adjustment to achieve realistic run lengths. 

A subsequent paper by the same authors [71J deals 
primarily with myosin VI, but it also presents a re- 
flned model for myosin V. This includes a simpler kinetic 
model, described with a 4-state cycle. The second impor- 
tant change is that the lever arm elasticity is modelled in 
a more realistic way, as an anisotropic elastic beam with 
persistence length values of 150 nm and 400 nm in two 
different planes. The elastic energy in a state where the 
two heads are separated by two actin subunits is reduced 
to about 10 kgT, as compared with 40 ksT in the origi- 
nal model. The calculated force- velocity relation remains 
similar as in the older model for forward and below-stall 
loads. For super-stall forces, the new model predicts 
nearly negligible backward velocities. 



S. Elastic lever arm model 

The elastic lever arm model |45j arose from a similar 
concept, but has many differences in the implementa- 
tion. The paper presents two alternative chemical cy- 
cles: the 4-state and the 5-state cycle. The 5-state cycle 
allows phosphate release from the lead head before the 
power-stroke, while the 4-state cycle does not. Note that 
model calculations and also later experimental evidence 
strongly favour the 5-state scenario, so we will focus our 
discussion of the model on it. The model uses a cycle 
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FIG. 7: Substeps in the elastic lever arm model. Because they 
involve changes between configurations bound with a single 
and with both heads, their size depends on load. 



of 4 attached and one detached state (plus the detached 
ADP state, which is not part of the regular cycle). Some 
short lived states (e.g., bound with ATP) are omitted, 
as are some unprobable ones (e.g., unbound head in the 
empty state). 

The lever arm is modelled as an isotropic elastic beam 
with a bending modulus of 1500pN/nm, which corre- 
sponds to a persistence length of 360 nm, or a spring 
constant of 0.25pN/nm measured at the lever arm tip. 
Both lever arms are connected with a fully flexible joint 
(the validity of this assumption was recently proven by 
Dunn and Spudich [IS]). 

Based on the above value for the elasticity of the lever 
arm we estimated that the lead head cannot commit its 
power-stroke before the trail head detaches, meaning that 
the telemark state does not occur. However, we cannot 
fully exclude that nonlinearities in the elastic properties 
of the lever arm still sometimes allow this configuration. 

The fact that ADP release takes place in the trail head 
while the lead head is attached to actin provides a direct 
explanation for the flat force-velocity curve below IpN 
load. Namely, although the ADP release is connected 
with a small conformational change in the trail head, 
it causes virtually no load movement. Conversely, the 
load therefore has no influence on the ADP release rate, 
which in turn determines the velocity. For higher loads, 
close to stall, other transitions (power stroke or lead head 
binding) become rate limiting and the velocity then drops 
abruptly. 

The model reproduces well the force-velocity relation 
up to the stall load. For super-stall forces, the model 
predicted backward shppage without ATP hydrolysis. It 
occurs when both heads have ADP on them. Very similar 
behaviour was later discovered by Gebhardt and cowork- 
ers |15j . However, the backward sliding speed in the 
model is too low. Another discrepancy is that experi- 



mentalists observe backward slippage even in the absence 
of any nucleotide, while the model needs at least ADP in 
order to move backwards. This is due to the assumption 
that the detachment rate of a head with ADP or with 
no nucleotide is the same for the lead and the trail head, 
and that it is independent of strain. In order to quantita- 
tively fit the newer data, a strain-dependent detachment 
rate in empty and ADP states will therefore be necessary. 

The model shows a short substep preceding the main 
step, which is connected with the trail head detachment. 
It also shows a sub-step following the main step, which 
corresponds to the binding of the lead head. The sec- 
ond substep is in good agreement with the observation 
of Cappello et al. [5^ . The first substep is not fully com- 
patible with their experimental result, because it shows 
no dependence on the ATP concentration. In total, one 
can still say that the substeps from that model are much 
closer to the experimental results of Cappello et al. [52] 
than to those of Uemura et al. [16^. In both interme- 
diate states before and after the main steps the motor 
is bound with a single head and therefore has a much 
higher compliance. Therefore, both substep sizes depend 
on the load. The first substep decreases with load (Fig. 
l7|, while the second increases. Uemura et al. [TB], on the 
other hand, see load-independent substep sizes. 

Another major advantage of the mechano-chemical ap- 
proach is that dissociation events that limit the run 
length are an inevitable and integral part of the model 
and do not need any additional assumptions or parame- 
ters. However, they do depend on several of the kinetic 
constants whose values are not yet well known to date. 
The model shows that there are generally 3 classes of 
dissociation pathways. 

Note that different models in the literature use quite 
different terminology for dissociation pathways. They are 
not entirely one-to-one translatable, mainly because the 
models of Baker j30] and Wu |48] include weakly bound 
lead head states with ADP.Pi and ADP on them, respec- 
tively. Other models do not distinguish between weakly 
and strongly bound states. So Table [Til presents the best 
effort of a "dictionary" , showing equivalent dissociation 
(termination) pathways in four models that include dis- 
sociation events. 

With the parameters used in the original paper, the 
model reproduces realistic run lengths for zero load. It 
also shows a run length that falls with ATP for concen- 
trations above 50 /^M. The run length decreases with an 
increasing ADP concentration - however, the findings of 
Baker et al. [30 are only partially reproduced. A more 
serious discrepancy is that the dissociation rate becomes 
very high for strong forward loads. Newer experiments 
have shown that myosin V stays attached to actin for a 
few seconds even under lOpN forward load |15j . There- 
fore, dissociation pathways 1 and 2 are likely overesti- 
mated in Ref. |l5] . 

The model was later expanded to take into account 
torsional fluctuations in the actin structure, which lead 
to an improved agreement between the calculated and the 
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TABLE II: Dissociation pathway "dictionary". The first column denotes the state of the dimer from which dissociation occurs. 
For example (A.M.T,A.M.D.Pi) means that the trail head is bound to actin with ATP on it, while the lead head is bound to 
actin with ADP and Pi. The other columns show the names of these dissociation pathways in different models. Evidence exists 
that "Path 3" is the predominant dissociation path. 



Pre-dissociation state 
(trail, lead) 


Vilfan |45] 


Baker et al. [30] 


Wu et al. [48] 


Skau et al. 47 


(A.M.T, A.M.D.Pi*) 


Path 1 


term3 
(excluded) 


terml 


- 


(A.M.T, M.D.Pi) 


Path 2 


- 


termmation 
(sole pathway) 


(A.M.D, M.D.Pi) 
(A.M.D, M.D*) 


Path 3 


terml 


term2 


- 


(A.M, A.M.D.Pi) 


- 


term2 


- 


- 



* weakly bound A.M.D in |i8] 



measured step size distribution 91j. It was also used to 
calculate the behaviour of myosin V when it encounters 
an Arp2/3 mediated filament branch j^l]. The predicted 
branching probability is in good agreement with data by 
Ali et al. ESI. 



V. CONCLUSIONS 

Our comparison has shown that there is no myosin 
V model so far that would fit all available experimen- 
tal data. A particular problem is to explain the high 
processivity over a wide range of loads. Additionally, 
our experimental knowledge about processivity is still 
somewhat incomplete. Data on load dependence of the 
run length for different ATP and ADP concentrations 
would be immensely useful. However, a glance at Ta- 
ble [n] possibly reveals the lines of a future consensus on 
the dissociation mechanism. Baker et al. |30j conclude 
that "terml" is the predominant dissociation pathway. 
In the elastic lever arm model 45J, three pathways are 
present, but only Path 3 shows a dependence that allows 
a high degree of processivity for forward loads. Wu et 
al. HH] consider two pathways, but use a higher disso- 
ciation rate for "term2". These terms all denote more 
or less the same dissociation pathway, namely dissocia- 



tion from the waiting, (A.M.ADP,A.M.ADP) state. The 
model by Skau et al. |47j . however, uses a completely 
different path. The consequence is that the run length 
strongly decreases with load, so that the model motor 
loses processivity when the load approaches IpN. 

With regard to sub-steps, there is still some discrep- 
ancy among experimental results [T5J [55] . This contro- 
versy will have to be sorted out. In addition, more data 
on load- and nucleotide dependence of the substep size 
and duration will provide valuable input for future theo- 
retical models. 

Finally, let us remark that all theoretical models de- 
scribed here are largely phenomenological. Connecting 
them in a quantitative manner with our knowledge of the 
conformational changes that take place inside the myosin 
head, such as opening/closing of switch 1, switch 2 and 
the cleft [27| [enilMl 113 [Ml EH [98], remains a major 
challenge for the future. 
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